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We compute the super-Hubble evolution of non-Gaussianity of primordial curvature perturbations 
in two-field inflation models by employing two formalisms: 5N and covariant formalisms. Although 
the two formalisms treat the evolution of fluctuations radically different, we show that the formulas 
of the /jvL parameter agree quantitatively with each other within 1% accuracy. We analytically find 
that the amplitude of /jvi decays no faster than a^'^ as the infiationary trajectory reaches to the 
adiabatic limit for generic potentials. 
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I. INTRODUCTION 

Cosmological perturbation theory is a powerful tool 
to analyze cosmological observations. It has been so- 
phisticated over years by confrontation with more and 
more precise data of the cosmic microwave background 
anisotropics [l|. 

In the single-field inflationary scenario the primordial 
curvature perturbations are generated during inflation 
and become the seed of the inhomogeneities seen today. 
These fluctuations are observed to follow nearly Gaussian 
statistics. It has been shown that all single-field infiation 
models with the canonical kinetic term produce almost 
Gaussian signals if the spectrum is nearly scale invariant 
[3-[5|[53j. The primordial curvature perturbations freeze 
in at the Hubble horizon exit, and the amplitude is con- 
served until the subsequent horizon re-entry regardless of 
the post-inflationary events This fact makes single- 
field inflationary scenarios very powerful and predictable; 
a flip side of this is that no one can know what happened 
after inflation and before radiation domination by ana- 
lyzing large-scale inhomogeneities on the sky. 

However, if curvature fluctuations are generated by 
some other fleld than the inflaton, they evolve until the 
inflationary trajectory reaches to the adiabatic limit (i.e., 
all entropy modes decay) and can be non-Gaussian 0- 
Q . Non-Gaussianities in multi-fleld infiation models have 
been studied extensively (see, e.g., [l0l - [33| ). In there, 
the key mechanism on producing non-Gaussian adiabatic 
perturbations always involves the existence of entropy 
perturbations on large scales. 

The 5N formalism [l^, 113, HI] is quite useful for 
the studies of non-Gaussianity since the curvature per- 
turbations can be obtained by computing derivatives of 
the local volume expansion (e-folding number N) with 
respect to the fields from initially flat hypersurface to 



finally uniform density hypersurface, where only back- 
ground equations need to be solved. 

On the other hand, in conventional perturbative ap- 
proaches to the spatial curvature, the perturbed Ein- 
stein equations and the perturbed Klein-Gordon equa- 
tions should be coupled and solved together with back- 
ground equations. Based on Hawkin g's covariant ap- 
proach to cosmological perturbations [36|, a nonlinear 
and covariant formalism for multi-scalar fields was de- 
veloped by (37l - l39j (see also [H, [l^ for an equivalent 
approach) to solve the evolution of curvature perturba- 
tions, and here we shall apply them to non-Gaussianities 
on large scales in the two-field infiation. Both formalisms 
are fully nonlinear and geometrically transparent, but 
treatments of the evolution of perturbations are radically 
different. Therefore, it would merit further study of non- 
Gaussianities by directly comparing the two approaches. 

In the 5N formalism, the non-Gaussianity is given by 
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where Jnl = /^l'''°" + /^T'^°'- Here iV/ = d/dip^ J Hdt 
and Njj = /dip^ d(f^ J Hdt are from background tra- 
jectories, while transfer functions T^'^^'*'^' represent super- 
horizon generation of curvature perturbations C from en- 
tropy ones. The subscript * denotes a value evaluated at 
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the Hubblc-horizon exit, and wc have used the standard 
slow roll parameters, e = —H/H"^ ,riij = V^i^j/V. 

In this note, we present a detailed numerical compar- 
ison between the two formalisms in case of two-field in- 
flation models. We show that Eqs. (IT|), © and ^ agree 
quantitatively within 1% accuracy. For illustrative pur- 
poses, we use canonical scalars with potentials used in the 
literature as examples, but the generalization is straight- 
forward and the same conclusion holds. A similar com- 
parison has been done by [13, |4l| for ekpyrotic models, 
but they have ignored the contribution from Eq. ([3]) re- 
sulting in the less accurate quantitative agreement. 

We also analyze the fate of /nl when the inflation- 
ary trajectory reaches to the adiabatic limit for generic 
potentials. In this case, we find that the amplitude of 
/nl converges to its asymptotic value no faster than 
contrary to the recent claim (22j . 

We employ the natural unit ?i = c = 1 and SttG = 



PI 



1 unless otherwise stated. 



II. COVARIANT PERTURBATIVE APPROACH 

In this section we review the evolution of primordial 
fluctuations on super-Hubble scales in the covariant per- 
turbativc approach. With the adiabatic and entropy field 
decomposition j37| (nonlinear generalization of [43|). one 
can clearly see how entropy modes convert into adiabatic 
curvature modes. 

For simplicity, we consider two light, canonical scalar 
fields 0, X in the Friedmann-Lemaitre-Robcrtson- Walker 
(FLRW) background: 
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Entropy perturbations up to the second order on large 
scales follow 



5s + 'iH6s + {Vss + 'i0'^)5s : 
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where Vsss = -(sine)^^^^^ -I- 'ico%B{smdYY^^^^ - 
3(cos6')2sin6'V^0^x + [coBdfV^^-^^. 

It is convenient to separate initial random Gaussian 
variables from their time evolutions as 

m = c* + 5sj:f\t) + 5siTf\t), (10) 

5s(t) = 5s^TfX^) + ^^lTfK^)^ (11) 

where and (5s* are evaluated at the Hubble exit. 

One can then write the solution of Eq. (|5]) with the 
transfer functions at each order: 



r(i)(t) 



dt' 



-2H 



-er, 



(1) 



i5s ' 



(12) 



-2H ■ 



(2) 



(7 



Ss 



(13) 



The super-Hubble evolution of the power spectrum of 
C is given by 
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We decompose their perturbations up to second order as 



5a = cos 05(j) + sin 95x + 
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where tan0 = x/cj) and a — + X^- Here dots denote 
derivatives with respect to the physical time. 

The adiabatic curvature perturbation is given by 



where ~ HI/Att'^. 

The super-Hubble evolution of the local-form bispcc- 
trum of C is given by Eqs. © and It Jnl becomes 
large after the Hubble exit, it is important to understand 
the behavior of f^f^^^'^'^ in detail. In fact, as we shall 
see, each term in T^^'^ [Eq. ^] contributes to 
equally, and thus one should not estimate it by select- 
ing a seemingly dominant term. Otherwise, it leads to 
a wrong conclusion that /jv^ is large since each of them 
gives a large positive or negative contribution if the en- 
tropic transfer occurs. 



C = --Sa + ^e5s5a -:^(-) 5a\ (7) 
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which is evaluated on a flat hypersurface. On super- 
Hubble scales, C evolves as 



III. TWO-FIELD INFLATION AND fNL: 
NUMERICAL ANALYSIS 



We employ a quadratic potential, such as 



C(i) w ~'^65s + ^{Vss + 4^2)^5^ - ^VJs5s, (8) 



mx) = ^0' + ^x'. 



(15) 



where = cos^F^ + sinOV^^ and Vss = (sin6')^V|00 — 
sin26'V;0x 4- (cos6')2t/^^. 



The inflationary trajectory of this model is characterized 
by a single parameter, a mass ratio R = m^/m^ [see 
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FIG. 1; A trajectory of double inflation with R = 1/9. 

Fig.[T]. We assume < m^. Also, we investigate other 
potentials in Appendix VK\ 

Non-Gaussianitics in two-field inflation with the po- 
tential (|15D have been quantitatively studied the most 
thoroughly in the literature. Vernizzi and Wands 
studied i? = 1/9 and Yokoyama et dX. ^ R = 1/20. 
Both of them used the 5N formalism and concluded that 
their results agree qualitatively with [l^ using the long- 
wavelength approach that is equivalent to the covariant 
approach. It was also shown that with various values of 
R the 5N formalism agrees numerically very well with 
the covariant lisi or the long- wavelength approach [l^ . 
Mulryne et al. j44l |45| showed that the momentum trans- 
port approach agrees with the 5N approach. 

In this section, we show that the covariant perturbative 
approach described above agrees quantitatively with the 
5N formalism for /jvl within 1% accuracy. 



A. Initial conditions and background trajectories 



In the background the system obeys 
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Inflationary trajectories are attractor solutions of coupled 
Eqs. (|16p and ([T7| : whatever the field values are chosen as 
initial conditions, the solution flow of the system quickly 
approaches to the attractor. We set the initial condi- 
tion on the attractor by solving the equations once and 
restarting from the attractor solutions. 

The initial conditions for perturbations are given by 
the linear perturbation theory [i^. A solution of the 
Mukhanov- Sasaki equation in the flat gauge is given by 
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and r\}-'^(x) 



is a Hankel 



where = 9/4 - tv? / H"^ 

function of the first kind. Here, hI^'{x <C 1) 
—iT{i'){x/2)~'^/7r is used in the last approximation. This 
approximation corresponds to the super-horizon limit, 
k ^ aH. Note that t = J dt/a{t) is the conformal time 
and lies in — oo < t < 0. Note also that the relation 
— fcr = k/{aH) is used above. 

We assume light scalar fields, then v — )■ 3/2, for which 
the spectrum is independent of k. We also assume that 
entropy modes must be light at the Hubble-horizon exit. 
Therefore, our initial conditions for fiuctuations are given 

by 
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Al (<5^i(ki)<5^;^(k2)5^f(k3)) =0, 



where S(pl is evaluated in the flat gauge. The last equa- 
tion implies S(p 
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0, however Sa^'' ^ and (5sl^^ ^ 



due to Eqs. jS]) and 

We assume slow-roll to give initial conditions for 
Eqs. (© and ®, and then Eqs. © and © reads [I^ 
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-M, (20) 



where e = a'^/{2H^), rj^^ = (cosfl)^//^^ + sm29r]^^ + 
(sin6')2?7^^, ri^^ (sin 9)'^ri^^ - sin26'?70^ + {cos 9)'^r]^^ 
and r]„s = sin20(77x^ - rj^^)/2 + (cos^ 9 - sin^ ^)?Mx- 

For the super-Hubble evolution, we calculate in two 
ways: (i) covariant formalism - we solve Eqs. ([U and (|9]) 
coupled with Eqs ((TE)) and p7|) from t = to some time 
after inflation and (ii) 5N formalism - we solve nine ad- 
jacent background trajectories with nine different initial 
field values, e.g., (^* + dip^x*), and then compute the 
first and second derivatives of the efolding number with 
respect to horizon crossing field values, such as 







N{<p + d^,x)~N{^-d^,x) 
2d(t) 
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where dip^ = 10~^ is chosen. We have checked that the 
results converged with 10"'^. For smaller values of dip^ , 
numerical noise tends to increase. Note that if we im- 
prove the treatment of derivatives by using the higher 
order of the difference than Eq. (HI]), numerical noise 
would decrease but it requires more trajectories in the 
calculation. 
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FIG. 2: Evolution of the power spectrum of which is nor- 
maUzed by = HII{A-k'^). Here A^(5Ar) indicates the 5N 



formahsm, while A? (trans) indicates the covariant formalism. 
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FIG. 4: Evolution of C^'^K The red line shows Eq. {THJ. The 
gray dash-dotted, blue dashed and pink dotted lines show the 
first, second and third terms in Eq. (|13|l . respectively. 
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FIG. 3: Evolution of non-Gaussianity, /jvl- Here fNL{SN) 
indicates the 5N formalism, while /jvi (trans) indicates the 
covariant formalism. 



For the SN formalism, we have to set spacetime folia- 
tion as initially flat and finally uniform-Hubble (comov- 
ing) hypersurfaccs. The curvature perturbations and the 
power spectrum are given by [3^ Issj 



C = SN = N^Scb, + N^Sx* + 



+N^^S(PJx* + ^N^xSr., (22) 

(23) 



Note that Eqs. ^ and ([9]) are gauge- invariant [s^l, thus 
we set the same foliation for a comparison purpose. 



double inflation with m^/m^ = 1/9 



We set the same initial condition as [l7[. Initial field 
values are = = 13. We assume that inflation ends 
at Ee = 1 and the number of efolds at the end of inflation 
from the Hubble-horizon exit is =60. We find that 
(0* = 12.9281, X* = 8.27202). The perturbed equations 
are integrated from this time. 

Figure [1] shows the trajectory in the field space, where 
one can find a turn. At the turn, entropy perturbations 
source C through Eq. Since the entropy mode be- 
comes heavy and decays sufficiently after the turn, it no 
longer sources (. As shown in Fig. [5] and Fig. |31 both 
amplitude and non-Gaussianity of C are generated dur- 
ing the turn. Both figures show agreement of two for- 
malisms with < 1% accuracy. The hight of the peak is 
^ 0.22 which agrees precisely in both the 6N formalism 
and the one calculated from Eqs. ([2]) and ([3]). Note that 
/nl computed here corresponds to in [l3|. Since the 
evolution of in [l3l is governed only by Nj and T^^^ 
for each formalism [i^ . the agreement is the same order 
as the power spectrum A^. 

How can we understand the feature of non-Gaussianity 
appeared during the turn? The cause of the tiny net 
effect is the precise cancellation between terms in Eq. ([8|) . 
As mentioned in Sec.|lll the shape of /nl is attributed to 
that of the second order transfer function [Eq. ([T5|) ] . 
As shown in Fig.|4l the second term has the sharpest rise; 
and then negative contributions from the first and the 
third terms catch up with it. The value after the turn is 
given by /wL - 0.007, which consists of /^T'^" ^ 0.005 



and /]jf°r"™" ^ 0.002. The value before the turn comes 
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FIG. 5; A trajectory of double inflation with R — 1/20. 



FIG. 8: Evolution of C^^^ 




FIG. 6: Evolution of the power spectrum of which is nor- 
malized by = H^/{4n^). 
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FIG. 7: Evolution of non-Gaussianity, /nl- The ratio is di- 
vergent (ill-defined) when /atl (trans) crosses zero. 



2. double inflation with m^jm^ — 1/20 

We set a similar initial condition to [l^]. Initial field 
values are set to (pi = Xi = 10, and we find ((/)* = 
9.99920, X* = 9.68230) by setting = 50. 

The similar feature as the case of i? = 1/9 can be 
found in Figs. [51 [HI and[71 There are a few sub-turns in 
a turn. Each sub-turn creates a series of peak and dip. 
Since the mass ratio is larger, the turn is sharper and 
leaves features on the power spectrum and bispectrum. 
Again, the SN and covariant formalisms agree with 
1% accuracy. The peak value is Jnl ^ 7.5. The value 
after the turn is /tvl ^ 0.001. The value before the turn 
is Inl ~ 0.002. One can clearly see the very delicate 
cancellation among the entropic transfer terms in Fig. [S] 

From the above two cases, one can gain some insights of 
the mechanism to generate non-Gaussianity at the turn 
in the field space. The sharper the turn is, the more 
non-Gaussianity is peaked. A peak is followed by a dip 
due to different timing contributions from each term in 
Eq. ([T3|: the net effect erases almost all the traces due 
to the attractor behavior after the turn. The residual 
non-Gaussianity can remain if the turn is slow enough to 
allow the background quantities, such as the expansion 
rate, change their values during the turn before entropy 
modes decay. 

We have checked that the accurate agreement also 
holds for other parameters and other potentials, e.g., 
V = mlx^e-^'''^ /2 used in [H Ell and F = 02%^ + 
bo — b2<fP + hicjy' used in [l^ (see Appendix E|. 



IV. ANALYTIC ESTIMATE OF DAMPING /a 



In the previous section we have shown that how the 
(2) 

cancellation of happens numerically when the en- 
tropy mode damps, i.e., the inflationary trajectory be- 
comes purely adiabatic since the entropy direction closes. 
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In this section, we analyze an envelope behavior of /nl 
when the trajectory approaches the adiabatic limit by 
solving the perturbed equations while leaving the form 
of the potential arbitrary. In this case, there are two 
regimes. Initially, entropy modes are light and source 
adiabatic modes. Eventually, entropy modes get heavier 
and damp away. We solve the super-Hubble evolution of 
C in the latter regime by dividing cases [A for ovcrdampcd 
(light) entropy modes, B and C for underdamped (heavy) 
entropy modes] , and we will find how fast Jnl approaches 
to the asymptotic value that can be large in principle. 

Physically, case A describes the slow-roll regime of 
two-field inflation, where curvature perturbations have 
already been sourced by the time ti > but entropy 
modes are still light. An example of this is double in- 
flation with R close to 1. Cases B and C describe the 
slow-roll regime, where entropy modes get heavy after 
sourcing Examples of these arc double inflation with 

= 1/4 for B and R = 1/20 for C. These three cases 
cover almost all two-field models when entropy modes 
decay while separating cases in actual models is subtle. 



A. Small r]ss and small 9 

When entropy modes are light {rjss < 3/4), one may 
use slow-roll equations for entropy modes to first order: 



out with the time evolution of (5s^^\ but we did not catch 
the feature in this estimation; we just obtained that the 
second order term damps away eventually if rjss is sizable. 
In this case, /nl C^^VCi damps as a"^''"" toward its 
asymptotic value. 



B. Large -qss and small 9 

In this case one cannot use the slow roll relations for 
entropy modes but can do along the trajectory, i.e., in 
CT-direction. 

From Eq. ([TS)) one has on the super-Hubble scale 
' ' V2H r(3/2)" [aHj 



30^ \ 



(28) 



where v is pure imaginary if r]ss > 3/4 and/or > 
3/4. Here, we shall assume large rjss > 3/4 while keeping 
other slow roll parameters small. This assumption leads 
to slow roll and slow turn, 9 ~ —{l + rjss/'i)i]asH. In this 
case, we have Ss^^^^ oc a~^^^. 

The super-Hubble evolution of ^ is given by 



,5s 



(24) 



C 



2H ■ 
'—0Ss 
a 



^(V^,,+4^2_ 9^2)^^2^ (29) 



where bs have been used in deriving the initial conditions 
C*, ,5s* in the previous section. We integrate these equa- 
tions further by assuming some entropic transfer hap- 
pened by the time ti > t^. 

The solution of Eq. ([21)) is given by 



(5s = 6s\ exp — / dtHr]s 
Jti 



const. 



(25) 



where we have used a slow roll relation rjss = 2Her]s 
27/772^ mill. 

The super-Hubble evolution of C is then given by 
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e e 
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where we have used slow roll relations: 'iHa ~ — Kr, 
(5s ~ —r/ssHSs and 9 j H ^ —rja-s — constant. Using 
Eq. ((25)) . one can integrate C as 



rjss V e V"i 



^-^s^f-V^T) 



where ^1 is the asymptotic value that is different from 
Eq. ()19p in general. The second order term should cancel 



where we have used 3H& ~ — and 5s 
Since rj^s ^ 2Herias + Hi]„s{v<ya~'nss)-Qo 
for r^ss ^ 1; we get 



- / dtHr]ss 











-{3/2)H5s. 
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= Vasl exp 

Using this and Eq. (|25)) . one can integrate ( as 
. , r?..i(77,,/3 + l) [2 /a\-''"-3/2 



Vs 



2e 



ai 



!?Li(W3 + i)^^^2 (± 

(77,, +3/2)e ^\ai 
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where we have assumed constancy of rjss- For ij^s > 3/4, 
the value of /nl decays as since the second line of 
Eq. (|3T)) dominates [see fig. [9]. 

Meyers and Sivanandam [22| have recently claimed 
that /nl decays exponentially when the entropy mode 
damps by using the 6N formalism (see Eq. (72) in 

Jnl ~ 0(e*) + 0{1) X 7j,s exp [-2 [ dtHC„ijJ , (32) 

where C,, is a number greater than 1 . This formula would 
correspond to the last term of Eq. (|M)) if ijss ^ const. ^ 
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FIG. 9: Damping of non-Gaussianity, log^Q I/jvl — /at^l™^ 
in the case oi R = 1/4, where f'^l"''^ = 0.00855163. The 
envelope scales as a^"^. The values of rjss and rjas in a linear 
scale are also shown for reference. 

3/2. This exponential behavior is valid only in the last 
tail of fNL that approaches to its asymptotic value and 
shows up after the oscillatory feature of entropy modes 
disappears. 

C. Large rjsa and large 9 

If the angle changes rapidly, one cannot use the slow 
roll relations for entropy modes or slow turn relations but 
still can do in cr-dircction. Instead, we use the scaling 
of ^ = (— 0^;^ + xV^0)/i7 ^ a^'^/^ since the heavy field 
decays as a~^/^ (if (j) is heavy, <p ^ a~'^/^, V^t, a~'^l'^ 
or faster depending on models). In this case, we can 
integrate Eq. (P^)) as 

where —9/H ^ Cela/ai)^^^^ have been used. Thus, we 
obtain Jnl - Vss{a/ai)~^ if Ce ^ 0{1) [see fig. [TO]. 



V. CONCLUSION 

We have re-examined the super-Hubble evolution of 
the primordial non-Gaussianity in two-field infiation by 
taking two approaches: the SN and the covariant pcrtur- 
bative formalisms. The results agree within 1% accuracy 
in the cases of two-field inflation with the quadratic po- 
tential ((15]) and with the other potentials in Appendix [XI 
The SN formalism and the long wavelength formalism, 
which is equivalent to the covariant one, were also com- 
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FIG. 10: Damping of non-Gaussianity, logj^Q |/jvl — /^fl^™''! 
in the case of i? = 1/20, where f^l"'"^ = 0.00854551. The 
envelope scales as a^'^. 

pared analytically and numerically in p^ . and our re- 
sults agree with theirs. Our treatment has improved the 
agreement in the previous literature. 

The peak feature appears on /jvi at the turn in the 
field space, which can be understood as the precise can- 
cellation between terms in Eq. ([TO]) . Note, however, that 
we have made the slow-roll approximation for both ap- 
proaches at the horizon exit even though integration after 
the horizon exit has been carried out without assuming 
slow roll. 

In Figs. [3] and [5] we have shown that the net effect of 
the second order entropic transfer is relatively small due 
to the cancellation between terms in Eq. ([TO]) with the 
quadratic potential, which also holds for the other po- 
tentials studied in Appendix [^ if the evolution is tracked 
until entropy modes decay. Thus, it is important to in- 
clude all terms in the second order perturbation theory 
for correctly estimating /nl- If /tvl is detected by ob- 
servations, the two-field model with Eq. ([TO]) is ruled out. 

Our method to calculate the second order ( can thus 
be used for cross-checking results with two formalisms, 
where careful evaluations of Eq. ([TO]) are necessary. In 
addition, inclusion of Eq. ([3]) is crucial to obtain the ac- 
curate agreement. 

We have analytically shown that if large Jnl is gen- 
erated during infiation, it decays as approaching to 
the adiabatic limit in a broad class of two-field models. 
In principle, persistent values of non-Gaussianity, f^]l™^ , 
may be observably large, but we did not meet such an 
example. 

Finally, let us mention the relative merits of the co- 
variant perturbative approach. This approach requires 
a single realization in the simulation; statistical proper- 
ties are taken into account in the formulation. It will 
significantly save a lot of computational time if integra- 
tion of the background equations is time consuming since 
the SN formalism requires many realizations, especially 
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FIG. 11: A trajectory of inflation model with Eq. (|A1|I . FIG. 12: Evolution of non-Gaussianity, /jvl for Eq. (|Aip . 



if one needs to deal with many scalar fields. It is straight- 
forward to gencrahze our formulation to N-fields. Also, 
the numerical noise is much less for the covariant formal- 
ism than the SN formalism. It still remains to be seen a 
computation of non-Gaussianity by the covariant pertur- 
bative approach coupled with a lattice simulation that is 
suitable for preheating. 
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Appendix A: Other examples: "large" 
non-Gaussianity 

In order to back up the results obtained, we further 
study two illuminating models which produce relatively 
"large" non-Gaussianities, I/tvlI ^ C'(l), in the sense 
that the values are much larger than those from single- 
field inflation and two-field infiation with Eq. (fTS)) . How- 
ever, since the upcomin g P lanck satellite is expected to 
constrain A/jvl ~ C(5) [5l|, these models do not predict 
Jnl large enough to be observed. Nevertheless, they are 
interesting examples. 

The models are 

Vi^,x) = ^X^e-^*^ (Al) 



V{(t>,x) = a2x' + &o - 62<^' + M^ (A2) 

which have been studied in [l^, 113 ^^'^ ' respec- 
tively. 

For Eq. (jAip . we set a similar initial condition to [l^, 
HI]. Initial field values are set to (0^ = 10"'^, Xi = 16), 
and we find = 0.00155720, x* = 15.4210) by setting 
Ne = 60. The model parameters are chosen to be A = 

0. 05, = 1 in numerical units to match with (27L l43|. 
Figs . 1 1 II and [T^ show that there are smooth, continuous 

turns generating non-Gaussianity of C until oscillations in 
X-direction start. Infiation ends at e = 1 and N = 60, 
which is right before the oscillation phase. The gener- 
ated /nl is negative and large, but transient. After the 
entropy mode gets heavy, the asymptotic value of /nl 
stays around /^^^'"'^ ^ —2.1 and a series of peak and 
dip shows up at each oscillation in x- Since each series 
of peak and dip at each oscillation cancels out precisely, 
/nl approaches to its asymptotic value as a~^. This 
fact was not pointed out previously in [i^ ; they stopped 
solving the evolution equations at the end of inflation. 

For Eq. (jA2p . we set a similar initial condition to [l6| . 
Initial field values are set to ~ 0.01, Xi = 18-1), 
and we find = 0.0100143, x* = 18.0239) by setting 
A^e = 85. The model parameters are chosen to be a2 = 

1, 62 = 7/20, 64 = 1/10,60 = ^2/(4^4) in numerical units 
to match with [31 . 

Figs . fT3l and fT4l show that there are a few turns gener- 
ating a few series of peak and dip in /nl- The first peak 
value is as large as fNL ~ 40. As in the case of Eq. (|A1[) . 
Jnl approaches to a persistent value of /^^]^™^ ^ 1.224 
as a~^. 

Tzavara and van Tent [l(i| have also compared the SN 
formalism and the long wavelength formalism, which is 
in principle equivalent to the covariant formalism after 
choosing a coordinate system and expanding to the low- 
est order in spatial gradients. In the case of Eq. (|A2p . 
they have reported asymptotic values of fj^^ ~ 1.19 for 
the long wavelength formalism and ^ 1.23 for the 
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5N formalism; they agree within 3.6%. Our formula for 

/nl [Eqs. © and ^] corresponds to Eq. (3.6) in 

In there, Vi2 corresponds to our ^y2e^T^^\ gint + <7iso to 

our —2t^T^\ W22 to our rjss, e — ??" to our r^o-o-, and —7;^ 
to our rjas in slow roll. 



FIG. 13: A trajectory ol inflation model with Eq. (|A2|> . 




FIG. 14: Evolution of non-Gaussianity, Jnl for Eq. (|A2|) . 
The ratio is divergent (ill-defined) when /jvl (trans) crosses 
zero. 
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In the following, we denote the local form of non- 
Gaussianity by fNL ■ Interestingly, the single-field consis- 
tency relation of /jvl holds also for non- Einstein gravity 
(see e.g., [12]). For the definition of /jvl, see [Slj . 
This solution can be also obtained from Eq. (|18|) by 
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